Effects of scattering area shape on spin conductance in a four-terminal spin-Hall setup 
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We study spin conductance in a ballistic and quasi-ballistic two dimensional electron system with 
Rasbha spin-orbit coupling. The system has a four-terminal geometry with round corners at the 
connection to the leads. It is found that by going from sharp corners to more round corners in 
the ballistic system the energy depended spin conductance goes from being relatively flat to a curve 
showing a series of minima and maxima. It is also found that when changing the size of the terminal 
area by modifying the roundness of the terminal corners the maxima and minima in the transverse 
spin conductance are shifted in energy. This shift is due increased (decreased) energy for smaller 
(larger) terminal area. These results were also found to be reasonably stable in quasi-ballistic 
systems. 

PACS numbers: 72.25.Dc, 71.70.Ej, 85.75.-d, 73.23.Ad, 72.20.Dp 



I. INTRODUCTION AND MOTIVATION 

Spin-orbit (SO) coupling in semiconductor nanostruc- 
tures has for some time been considered as one of the 
main candidates in controlling spin in semiconductor 
spintronic devicesA^-^i. One of the main benefits of the 
SO coupling is that it can be used to manipulate electron 
spins in semiconductors through gate voltages instead of 
external magnetic fields^ or magnetic doping-. SO cou- 
pling arises from relativistic effects and when combined 
with k ■ p - calculation it will lead to spin-orbit terms 
depending on the crystal structure, i.e. the Dresselhaus 
SO coupling^'^ in crystals that lack inversion symmetry. 
In addition, further spin-orbit contributions can occur at 
heterostructure interfaces. This contribution is the so- 
called Rashba SO coupling^^iS. Voltages applied to local 
gates can change the heterostructure confining potential, 
thus modifying the Rashba SO coupling^' More re- 
cently, yet a new type of spin-orbit interaction has been 
found in symmetric two-dimensional quantum structures 
with two subbands: the inter-subband induced spin-orbit 
couplingJiiiS,. 

In the field of spintronics the main goals is the cre- 
ation and detection of spin-currents. In electron doped 
semiconductor structures one of the main candidates for 
spin-current creation is the so-called spin-Hall eSecti^. 
Although it has been proved that in extended 2D sys- 
tems the spin-Hall effect vanishesi^"— , this results does 
not hold in finite size systems. The exception here 
are systems with the inter-subband induced spin-orbit 
coupling^ ^'1^, which gives rise to a non-zero spin Hall 
effeclj^, even for extended systems. The interplay of 
in-plane confinement and spin-orbit interaction leads to 
many interesting spin-related transport phenomena. 

The Rashba SO coupling has been in recent years ex- 
tensively studied both experimentall y ^^d^" — and theo- 
retically, either in two-terminal setups^^— , or in multi- 
terminal setups21i^^— . Previous numerical studies hint 
that the shape of the scattering area of the system plays 
a vital role in the behavior of the spin current fiowing 



through the transverse leads^i^. Also one expects that 
gate defined nanostructures will have confinement poten- 
tial that changes on a length scale much larger than the 
Fermi wavelength of the electrons"^^ . 

In this work we will present calculation of spin conduc- 
tance through transverse leads in a four terminal spin 
Hall setup where the shape of the scattering region is 
changed. Most lattice model calculation o^^'^"" — used 
abrupt edges in the way the leads where connected to the 
scattering region. This abruptness leads to substantial 
scattering (both of charge and spin) which can suppress 
spin-related phenomena that one is interested in. We 
use smooth connection of leads to the scattering area. 
Since most semiconductor heterostructures are defined 
by gates, one would expect relatively smooth confining 
potentials^. We study the spin conductance as a func- 
tion of scattering area shape and propose an explanation 
of the observed connection between the spin conductance 
and the change of size of the scattering region. We also 
observe that the values of spin conductance is relatively 
large which is due to the smooth transfer of electrons 
from the longitudinal leads to the transverse one. 

The paper is organized as follows. In section II we 
present the theory behind the calculations. Section 
HI defines the parameters used in the calculations and 
presents the results for the ballistic system in subsection 
A and for the quasi-ballistic system in subsection B. Fi- 
nally section IV contains conclusions and discussion. 



II. THEORY 

We are interested in a four terminal 2DEG in a semi- 
conductor heterostructure. A schematic of the scattering 
region of the system can be seen in Fig. [1] The system 
is described by the following effective mass Hamiltonian 



pI+pI 

2m* 



+ Vc{x,y) + Hn{x,y), (1) 



where m* is the effective electron mass. Vc{x,y) de- 
scribes a hard-wall confining potential shape of the four 
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terminal junction with round corners. The roundness 
of each corner can be controUed independently by vary- 
ing the radii of the circular corners, Rqi and i?c2- In 
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FIG. 1. Schematic picture of the four terminal system. Sharp 
corners are replaced by round corners with radii Rci and Rc2, 
see Fig. [1] The whole system has a length of Lx = 102a = 600 
nm and width of Ly = Lx- Leads 1 and 2 are the main leads 
into the system with a width of Wy = 41a = 241 nm, while 
leads 3 and 4 are the transverse leads into the system with a 
width of Wx ~ 18a = 106 nm. The Rashba SO coupling is 
turned on smoothly on a 10a = 58.8 nm thick strip framing 
the system (hatched area in the figure). Inside the scattering 
region is shown a circle of width D which we use to estimate 
the size of the terminal area. 



function. All the matrices are 2NxNy x 2NxNy matrices. 
In order to save computational power only the necessary 
Green's function matrix elements are calculated with the 
recursive Green's function method'^^. From the Green's 
functions both the spin conductance and spin-densities 
can be calculated. The z spin density for an bias of cVq 
is 

(S.(r)> = - / :^Tr,pi„ [a.G< (r, r; i?)] , (5) 

where is the Pauli matrix for the z direction. 



G<{r,r;E) = i 



is the lesser Green's function 
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(6) 

-2Im(Sj), and 



fj{E) = l/[exp(/3(£' — fij)) + 1] is the Fermi function in 
lead j which has the chemical potential — Ep + eVj . 
We define the spin-dependent conductance between the 
leads q and p as 



h 



(7) 



Using the Landauer Biittiker formula we write the linear 
response of the a ^ spin depended current through 
lead p as 



(8) 



a region centered on the scattering area we have a SO 
coupling which is described by the symmetrized Rashba 
Hamiltonian 

Hr{x, y) ^ {aix, y), {py(T^ - Px(Ty)} ■ (2) 

Here a(a:, y) is the Rashba SO coupling strength which 
is turned on smoothly at the edges of the scattering re- 
gion. Having Rashba SO coupling only in a finite area 
of the sample can be achieved by using metallic gate on 
top of the 2DEG'^''. The Hamiltonian is discretized on a 
square lattice using the finite difference method, resulting 
in an infinitely large matrix. The corresponding matrix 
equation for the retarded Green's function matrix of the 
system is 

{El - H)G'' = I. (3) 

The scattering region part of the Hamiltonian Eq. ([3]) can 
be made finite by treating the leads as self-energies in the 
standard way2L2§,. The leads have no SO coupling. The 
finite version of the matrix equation is 

(ii;i-H,-^S,)G^ = I. (4) 

Here the self-energy of lead j is S j , Hg is the Hamiltonian 
of the scattering region and Gg is the retarded Green's 



III. NUMERICAL SIMULATIONS 

The dimension of the scattering area is set within a 
Lx X Ly box shaped area, where = Ly = 600 nm. 
This area is discretized on a 102 x 102 point grid with 
constant mesh size a = 5.88 nm and the effective electron 
mass is set as m* — 0.0447me, which is a reasonable 
value for Gai_a;Ina;As alloys^. In the simulations all 
lengths where scaled in the mesh size a and energies in the 
tight-binding hopping term t = Y? I2m*a} = 24.6 meV 
which results from the finite difference discretization of 
the Hamiltonian in Eq. ([T|). Four leads connect to the 
system, two main leads with width Wy = 41a = 241 
nm (leads 1 and 2) and two transverse leads with width 
14^2: = 18a = 106 nm (leads 3 and 4). In the following, the 
energy will be scaled in the lowest transverse energy of 
the main leads, i?o = 0.145 meV= 5.89 • 10~^t. Similarly 
we also define E^ = 5.46£^o which is the lowest transverse 
energy of the transverse leads. The effect of Rashba SO 
coupling is turned on smoothly over a 10a = 58.8 nm 
strip at the edge of the scattering area, see the hatched 
area in Fig. [1] to its full strength a^. Between the two 
main leads we apply a bias eVo — IEq. This means that 
the chemical potentials at the leads measured from Ep 
will be eVi = eVb/2, eV2 = -eVb/2 and eVg = el^ = 0. 
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According to Eq. ([8]) the spin depended current {a =t, i) 
is then 

= (GS^ + G^f - G?r - G^f ) y t = 3,4. (9) 

We define the spin Hall conductance in transverse lead t 
as 

G.H. = :^^^ i = 3,4. (10) 

Note that in our scenario the total charge current through 
the transverse leads is always zero because of how we de- 
fine the bias over the leads. Throughout the paper the 
full strength of the Rasbha SO coupling will be set as 
ao = 10 nieV nni= 2atso, where ^so = 3.46 • IQ-H. This 
value of Rashba SO coupling is typical for Gai^^^In^jAs al- 
loys and is also in the regime where we begin to see strong 
and interesting features in spin conductance through the 
transverse lead. 



A. Ballistic system 

We will start by showing the spin Hall conductance 
through the transverse leads in a clean system. In Fig. [5] 
the spin Hall conductance through lead 3, GsHs, is plot- 
ted as a function of Fermi energy E-p for three types of 
shapes. Fig. [5] a) corresponds to shape with Rq — 29a, 
Fig. [2]b) to Rc — 26a. In both figures the the spin con- 
ductance for a sharp corner system (i?c = Oa) is plotted 
for comparison. Here we use Rq = Rci — Rc2 when all 
the corners have the same curvature. 

In Fig. [5] a) and b) we see that the spin conductance 
curves corresponding to the round corner system shows 
a series of minima and maxima. These extrema are most 
likely related to the resonance states in the scattering 
area which we will discuss below. 

The resonance states form when the leads get adia- 
batically larger which results in a downward shift in the 
transverse energy bands of the leads'*^, see Fig. [3] This 
causes new states to become available at the Fermi en- 
ergy which have very low longitudinal velocity {v^ ~ 0). 
Near the transverse leads the adiabatic approximation 
ceases to apply and the electron scatter in to different 
states, including the newly opened low longitudinal ve- 
locity states which give rise to conductance extrema in 
the transverse leads. 

The sharp corner GsHs curve, shown both in Fig. [5] 
a) and b), behaves differently. At high enough energies 
[Ep > 30i?o) the curve changes slowly with energy and 
has no significant resonance peaks, apart from isolated 
peaks due to divergence in the density of states at ener- 
gies corresponding to the band bottom in the longitudinal 
leads. The observed behaviour of the sharp corner GsHs 
is due to the strong scattering by the sharp corners. An 
electron coming in from the left, with a definite /c-value, 
gets scattered into all possible states, with the same en- 
ergy, when it enters the scattering region. For higher 



Rc=29a R^.=26a Rc= Oa 




Ep [Eq] 



FIG. 2. (Color online) Spin conductance through lead 3 
plotted as a function of Fermi energy. Figure a) correspond 
to the four terminal shape where Rc = 29a, b) correspond 
to Rc = 26a. In both figures the spin conductance through 
a Rc ~ Oa system (sharp corners) is plotted for comparison. 
The bottom of the energy bands of transverse leads (with 
energy scaled in Et) have been marked onto the figure. Note 
that the y scale is unique for each of the figures. 

energies these states are many, all corresponding to dif- 
ferent effective magnetic fields which tend to average out 
the spin-signal. Note that for low energies we still see 
peaks and minima since only a few fc-states are available 
at such low energies. 

Comparing the round corner curves in Fig. [5] a) and 
b) we see that the extrema at Ep = 56Eo in Fig. [2] a) 
has shifted to Ep = QOEq in Fig. [2]b). Other extrema 
in GsHg seem also to have been shifted by the change in 
corner radius Rq. To explore in more detail this shift 
of extrema in GsHa with change in curvature we make a 
surface plot of the spin Hall conductance through lead 3, 
GsHa J as a function of Fermi energy varied from OEq to 
100i?o and Rc varied from Oa (sharp corners) to 30a with 
interval of la. The result can be seen in Fig. |4]a). There 
we see that for low corner curvatures {Rc < 10a) the 
spin conductance becomes smeared out as was discussed 
above. 

For higher corner curvatures (Rc ^ 10a) we get a spin 
conductance curve which shows a series of minima and 
maxima. Extrema in GsHs also begin to shift in energy 
with increasing corner curvature. For smoother corners 
the transport gets more adiabatic and the extrema in 
GsHa will be more influenced by geometric resonances, 
i.e. changes in shape will be adiabatically translated into 
a shift in the energy bands, that will affect the trans- 
port. Changing the curvature of the scattering region ef- 
fectively changes its area. For smooth enough corners, the 
shape of the system is not affected by changing the cur- 
vature, only its size, which results in a universal shift of 
all bands in the scattering area. This does not apply to 
the sharp corners, where a small change in curvature can 
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results in great change in scattering properties. 

The effective diameter of the scattering region is D — 

{V2 - l)(i?ci + Rc2) + + W2, see Fig. □ The 

energy shift of a given GsHs extrema will be E(j cx D^^. 
The prefactor is not known, but this does not matter 
as we will show here below. We can pick a reference 




FIG. 3. (Color online) Schematic picture showing the effec- 
tive shift in the transverse Rasbha energy bands, see a), when 
the system gets wider, see b). In a) the left energy bands cor- 
responds to the system where the width is Wi and the right 
energy bands correspond to the system where the width is 
Wii. Here kx is the longitudinal k- vector of the electrons. 
The red lines correspond to spin up states and the blue lines 
to spin down states. At Fermi energies Ei and E2 we get addi- 
tional states with low longitudinal velocities (wx « 0), marked 
with white dots, in addition to the states with positive veloc- 
ities, marked with red and blue dots. In the terminal area 
the electrons from the original positive longitudinal velocity 
states can scatter into these newly open states. 



energy Eco, with reference size Dq = (v^ — l)(i?cio + 
i?C2o ) + V^A + same extrema in GsHs ■ By 

dividing Ec with this reference energy Eqq we get rid of 
the unknown prefactor in Eq and obtain 
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(V2 - l)(i?ci„ + Rc2„) + JW^ + Wl 



(V2 - l)(i?ci + i?c2) + JWl + Wl 



E, 



Co 7 



(11) 

which describes the shift of the energy levels in the scat- 
tering area as a function of the shape of the scattering 
area. 

In Fig. U a) we have plotted curves, shown with black 
dashed lines, with reference points at i?Co = -^Cio = 
^C2o = 30a on all the spin conductance maxima. As can 
be seen for Rc ^ 10a the shift of the maxima curves 
fits well to the shift that we expect for adiabatic change 
in size of the scattering area. For Rq < 10a scattering 
from the corners gets more dominating, which is to be 
expected since the system is out the adiabatic regime 
and the sensitivity to the corners takes over. 

A shape where we keep i?c2 = 30a fixed and only vary 
Rci was also examined, see Fig. |4]b). Here we also plot 
curves based on Eq. ()11|) . shown with black dashed lines, 
with reference points at i?ci = 30a on all the maxima 



curves. We notice that the for low i?ci corner curvatures, 
Rqi < 10a, we get somewhat more structure than for 
low Rq in the equal corner shape. This is due to that 
only two corners are contributing to the sharp corner 
scattering. We also notice that for i?ci ^ 10a the shift 
in spin conductance maxima follows the expected shift 
calcuated from Eq. (fTT]) . 




10 20 30 40 50 60 70 80 90 100 

Ep [Eq] 

FIG. 4. (Color online) Spin conductance through lead 3 
plotted as a function of radii of the corners and Fermi energy. 
In a) all corners are equally varied while in b) two corners 
are varied equally, the Rci corners, and two corners are kept 
fixed, Rc2 = 30a. The black dashed lines show the shift of 
GsH maxima in energy according to Eq. Hi}. The horizontal 
black lines mark where the spin conductance through lead 3 
in Fig. [2] a) and b) lie and the black circles show the (-Rc, ^^f) 
coordinates of the spin densities shown in Fig. [5] The vertical 
black lines mark the bottom of the transverse energy bands 
scaled in Et = 5A6Eo. 

To examine the effect of corner curvature better we also 
plot the density of the z spin component in the scattering 
region for three types of all equal corner shapes Rc = 
29a, see Fig. [5] a), Rq = 29a, see Fig. [5]b), and Rc = Oa, 
see Fig.[S]c). For all three types we set the bias as eVo = 
IEq and temperature at T = IK. 

In Fig. [5] a) we have stronger spin densities along the 
corner edges and within the transverse leads than in Fig. 
Ob). Also the spin densities averaged over the transverse 
leads in Fig. [S] b) is lower compared to equivalent aver- 
aging in Fig. [5] a). These densities correspond to spin 
conductance strength marked with black circles in Fig. |4] 
a) and as expected the more negative spin conductance 
yields more negative spin density averaged over the trans- 
verse leads. For comparison we also include spin density 
in the sharp corner system, see Fig.[5]c). We see that the 
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spin density in the sharp corner system is more smeared 
out than in the round corner systems, as could be ex- 
pected in hght of previous discussion about the sharp 
corner system. Note though that these spin densities 
seen in Fig. [5] are not large, only up to ^ 1.3 • 10^^ h/2 
per nm^. 
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1. Static impurities 

Phenomena such as crystal defects can introduce extra 
potential bumps in the otherwise uniform potential back- 
ground. These bumps can be described as static Gaus- 
sian shaped impurities in the scattering region which we 
add randomly to our scattering region via the term 



Ni 



Vi{x,y) = ^exp 



xf + {y^ - vf 

2Vl 



{12) 



which is added to the Hamiltonian in Eq. ([T}. Here Vj^ 
is the potential height of the potential bump, A^i is the 
number of impurities, {xi,yi) the center point of each 
impurity and Fj impurity range. We would expect the 
number of these impurities to be just a few percent of 
the number of donors. We choose 2% which corresponds 
to roughly Ni — 20 Gaussian impurities in an 600 x 600 
nm^ area and a delta donor density of w 3 • 10^^ cm~^. 
We are interested in seeing how strong these Gaussian 
shaped impurities must be to have significant effect on 
the spin hall conductance. To accomplish this four cases 
of impurity strength were studied: Vi^ = li?o, = 5i?o, 
Via — 10-Eo, and Vi^ = 50i?o, see Fig.O These values of 
impurity strength are roughly the same order of magni- 
tude as a screened and unscreened point charge 17 nm 
away from a 2DEG in Ga[Al]As heterostructure'^^. The 



FIG. 5. (Color online) Density of the spin z component in the 
ballistic system at _Bf = SB-Eo for three different geometries. 
The shape in a) has Rc ~ 29a, the shape in b) has Rc = 26a, 
and the shape in c) has Rc ~ Oa (sharp corners). These 
{Ep,Rc) values correspond to the spin Hall values marked 
with black circles in Fig. |4] a) . 

The results presented in this section suggest that po- 
larization of the spin current through the transverse leads 
can be effectively controlled by tuning the curvature of 
the corners in the scattering area. This could be realised 
by using e.g. finger gates^. 



B. Quasi-ballistic system 
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To test how robust the Gsh extrema are, we add im- 
purity effects to our model. Two different methods of in- 
cluding impurity effects were studied. The first method 
which we will consider involves adding a few randomly 
distributed Gaussian shaped impurities to the sample. 
The second method which we considered is the Ander- 
son impurity method which has already been extensively 
used for similar simulations^i"— i^i^. 



FIG. 6. (Color online) Scattering shape using Rc ~ 29a 
corner radii and with A'l — 20 static impurities added on 
random locations in the whole box of the scattering area. All 
impurities are set to have range Fi — 2a, and the impurity 
strength is set at Vig — IEq for impurities in a), Vig — 5Eo 
for impurities in b), Vig — WEq for impurities in c), and 
Vif, = 50i5o for impurities in d). 

effect these impurities have on the spin Hall conductance 
in leads 3 and 4 can be seen in Fig. [T] There we see that 
up to impurity strength Vig = IOEq the spin conductance 
is rather stable and that spin conductance through lead 
3 and 4 is symmetric around zero. 
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Vi=OE, 



Vi°=5Eo 







Vi=10Eo 
Vi=50Eo 




10 20 30 40 50 60 70 80 90 100 
Ep [Eq] 

FIG. 7. (Color online) The spin conductance through lead 
3 and 4 as a function of energy for the cases shown in Fig. 
[5] compared to the spin conductance of the clean system, see 
Fig. [2] a). Note that the Gsh scale is shifted incrementally for 
each Vifj curve from the Vi^ = 1-Bo curve. 



2. Anderson type impunties 

The Anderson impurity potential adds to each point 
in the scattering area a random value of potential en- 
ergy in the range [—U/2,U/2], where U describes the 
disorder strength. To study the effect of the Anderson 
type impurities on Gsh we focus on a extrema in GsHs 
at Ep = 56Eo for a Rq = 29a scattering shape, see Fig. 
[2] a) . In Fig. [8] the spin Hall conductance is plotted for 
impurity strength ranging from U = 0.05t to U = l.Ot. 
Each data point is averaged over 1000 impurity config- 
urations. For comparison we include the result for the 
ballistic system, seen also in Fig. [2] a). In Fig. |8]we see 
that the structure of the spin conductance is stable up 
to U ~ 0.2t = 3A.2Eq although its amplitude diminishes 
and the structure shifts with increasing disorder strength. 
For U = OM = 85.5i?o the spin conductance seems to be 
disappearing as has been shown for such large disorder 
values^. 



the ballistic regime and the quasi-ballistic regime were 
considered. It is found that above a certain curvature 
the spin conductance, Gsh, shows a series of minima and 
maxima, as a function of energy. These minima and max- 
ima are shifted in energy with increasing curvature. We 
propose a relation between this shift of the extrema and 
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FIG. 8. (Color online) The spin conductance through lead 3 
as a function of Fermi energy for a range of disorder strengths 
U . To compare we also plot he spin conductance of the bal- 
listic system which we mark with U — Ot. These results have 
been averaged over 1000 runs and the disorder strength is 
scaled in the tight binding hopping term t = 171i?o. 



the change in the size of the scattering area when the cur- 
vature is increased. This is due to the increased energy 
for smaller scattering areas and describes the observed 
shift of the extrema reasonably well. With this relation 
the polarization of the spin current through the trans- 
verse leads could effectively be controlled by tuning the 
curvature of the corners in the scattering area. 

In addition, we have also test the robustness of the spin 
conductance extrema by adding impurities to the scatter- 
ing area. Static and Anderson type impurities were con- 
sidered. The amplitude of the spin conductance extrema 
are suppressed with increasing impurity strength but the 
structure is shown to be stable under reasonable values of 
impurity strength, i.e. when the system is quasi-ballistic. 



IV. CONCLUSIONS AND DISCUSSION 

In this paper the effect of scattering area shape on 
spin conductance in a four-terminal 2DEG structure with 
Rashba SO coupling is studied numerically. The change 
in the scattering region was implemented by modifying 
the curvature of the corners between different leads. Both 
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